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Abstract 



Astrophysical jets have been studied with observations, theoretical models and numerical 
O ^, simulations for decades. Recently, supersonic magnetized jets have been formed in laboratory 

I ■ experiments of high-energy density plasmas. I will review these studies and discuss the 

O ' experimental setup that has been used to form millimeter-scale jets driven by strong toroidal 

magnetic fields in a MAGPIE generator. The physical conditions of these experiments are such 
^ ■ that they can be scaled to astrophysical scenarios. These laboratory jets provide insights on 

the underlying physics of magnetic tower jets and help constrain some models of astrophysical 
jets. In this context, we also discuss the connection between the laboratory jets and recent 
3D-MHD numerical simulations of Poynting flux dominated jets. The simulations allow us 
^ \ to investigate the effects of thermal energy losses and base rotation on the growth rate of 

kink mode perturbations, and to compare the evolution of PFD jets with a hydrodynamic 
counterpart of the same energy flux. 



1 Introduction 

Jets are observed in Young Stellar Objects, post-AGB stars, X-ray binaries, radio galaxies 
and other astrophysical objects. Models suggest that jets are launched and collimated in- 
side their "central engine" by a symbiosis of accretion, rotation and magnetic mechanisms 
(Pudritz et al. |2QQ7p . The engines cannot be directly observed however, because telescopes 
have insufficient resolution. Recently, laboratory astrophysics experiments have provided 
scale models of the launch and propagation of magnetized jets with dimensionless parame- 
ters relevant for astrophysical systems (Lebedev et al. 2005; Ciardi et al. 120091 Suzuki- Vidal 
et al. I2010p . These studies, when combined with numerical simulations (Ciardi et al. I2007( 
Huarte-Espinosa et al. I201ip . can help to resolve unanswered questions in jet physics. A 
fundamental distinction can be made between the physics of jet launch and that of jet propa- 
gation far from the engine. Since we cannot resolve the former observationally, it is important 
to identify distinct features of jets in the asymptotic propagation regime that can distinguish 
different engine paradigms. While both simulations and experiments now consistently reveal 
the promise, if not essentiality, of dynamically significant magnetic fields for jet launch, the 
correlation between the initial jet magnetic configuration and the stability of the flow far from 
the launching region is unclear. The importance of the magnetic field relative to the flows' 
kinetic energy divides jets into (i) Poynting flux dominated (PFD; Shibata & Uchida |T986j) , 
in which magnetic fields dominate the jet structure, (ii) magnetocentrifugal (Blandford & 
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Payne [T982j) , in which magnetic fields only dominate out to the Alfven radius. The observ- 
able differences between PFD and magnetocentrifugal jets are unclear, as are the effects that 
cooling and rotation have on PFD jets. 

2 Laboratory experiments 

Experiments and numerical modelling of astrophysically relevant supersonic plasma jets have 
been performed by a number of authors using high intensity lasers (see e.g. Farley et al. ll999| 
Shigemori et al. I2QQQ| Foster et al. I2QQ2( see also review in Blackman 2007). In this sec- 
tion however, we focus on magnetized jet experiments which have dimensionless numbers 
(Reynolds, magnetic Reynolds and Peclet) in reasonably appropriate regimes for astrophysics 
and have used conical wire arrays and conducting foils on pulsed power facilities (Lebedev et 
al. 2005; Ciardi et al. 2007,2009; Suzuki- Videl et al. 2010). 

2.1 Wire array, toroidal field 

Lebedev et al. (^2005^) applied a TW electrical pulse (1 MA, 250 ns) to an array of wires located 
inside a vacuum chamber. The set up consisted of a pair of concentric electrodes connected 
radially by tungsten wires of 13 /im in diameter. The current causes ablation of the wires 
which results in the formation of a background ambient plasma (Fig. la). This material is 
then pushed above the wires by Lorentz forces, and resistive diffusion keeps the current close 
to the wires. The current induces a toroidal magnetic field which at this stage is confined 
below the wires, around the central electrode. Then, after the complete ablation of the wires 
near the central electrode, the current switches to the plasma and creates a magnetic cavity 
(Fig. lb) with a jet at its center. The core of the jet is confined and accelerated upward by 
the pressure of the toroidal field. The return current flows along the walls of the magnetic 
cavity, which is in turn confined by both the thermal pressure and the inertia of the ambient 
plasma. Next, the magnetic cavity opens up, the jet becomes detached and it propagates 
away from the source at a velocity of order 200kms~^ (Fig- Ic). 
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(a) (b) (c) (d) (e) (f) 

Figure 1: Time sequence of soft X-rays images showing the formation of the background 
plasma (a), the expansion of the magnetic cavity (b), the launch of the jet (c), the development 
of instabilities in the central jet column (d) and the fragmentation of the jet (e,f). The wire 
array shows in the bottom part of the figures. Image taken from Suzuki- Vidal et al. ()2010p . 

Finally, instabilities, which resemble those of the kink mode (m =1), develop within the 
body of the jet (Fig. Id). The outflow is then fragmented into well collimated structures with 
characteristic axial non-uniformities (Fig. le,f; Lebedev et al. 12005^ . We note that a critical 
ingredient of these experiments is the significant thermal energy loss of both the jet and the 
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ambient plasmas. This is relevant because cooling plays a critical role in many astrophysical 
jet environments, e.g. YSOs. 

The experiment was reproduced a few years later by Ciardi et al. (12007 ^) using 3-D non- 
ideal MHD simulations. A numerical model was carefully designed to simulate the laboratory 
components (electrodes and wires) and all the plasma evolution phases. Good agreement was 
found between the simulations and the experiments (Ciardi et al. I2QQ7|) . The simulations 
showed that during the final unstable phase of jet propagation, the magnetic fields in the 
central jet adopted a twisted helical structure. This confirmed that the jets are affected by 
normal mode, m = 1, perturbations. 

2.2 Wire array, toroidal and axial fields 

A subsequent series of experiments were carried out by Suzuki- Vidal et al. (|2QlQp who studied 
the effect of introducing an axial magnetic field, 5^, into a radial wire array. Their experi- 
mental configuration was the same as that of Lebedev et al. (2005^ above), but they placed 
two solenoids below the wires, into the path of the current. Suzuki- Vidal et al. (120 lOp found 
that the added Bz affects the degree of compression of the on- axis pinch of the jet. As they 
increased the magnitude of B^^ the plasma column radially expanded after reaching a mini- 
mum radius. They also saw that jets are more stable when Bz is present than otherwise. We 
note that this stability effect can be explained analytically (with perturbation theory on a 
plasma column) and numerically (section 3) as a result of magnetic pressure form Bz acting 
against hoop stress compression from the toroidal field component. 

2.3 Conducting foil 

Ciardi et al. (|2QQ9 ) extended the experiments of Lebedev 
et al. ()2QQ5p by replacing the wire array (section 2.1) by a 
6 /im thick aluminum foil. Ciardi et al. ()2QQ9p found the 
following. The conducting foil experiment produced very 
similar results than the wire array one. However, while 
the wire array experiment results in the launch of one jet 
only, the conducting foil experiment results in a series of 
jets which are launched sequentially (Fig. 2). The plasma 
flux caused by ablation of the foil is much larger than 
that of the wires; the foil provides an increased mass as a 
function of radius. Thus the current gap produced by the 
magnetic cavity (section 2.1) is smaller in the foil exper- 
iment than in the wire array one, and can be refilled by 
the readily available plasma. Closure of the gap, which 
does not happen in the wire array case, allows the cur- 
rent to fiow once again across the base of the magnetic 
cavity, thus re-establishing the initial configuration. Once 
the magnetic pressure is large enough to break through 
this newly deposited mass, a new jet/bubble ejection cy- 
cle begins (Ciardi et al. 2009). This constituted the first 
laboratory study to address time-dependent episodic jet 
launch. 
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Figure 2: Filtered XUV emis- 
sion images of one of the episodic 
jet experiments. The clumpy jet 
at the center is contained inside 
three nested magnetic bubbles 
which were formed by previously 
launched jets (Ciardi et al. I2009p . 
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3 Simulations 

3.1 Model 

We use the Adaptive Mesh Refinement code AstroBEAR2.0 (Cunningham et al. 120091 Carroh- 
Nellenback et al. i201ip to solve the equations of radiative-MHD in 3D. The grid represents 
160x160x400 AU divided into 64x64x80 cells plus 2 adaptive refinement levels. We use 
periodic boundary conditions at the four vertical faces of the domain, extrapolation conditions 
at the top face and a combination of refiective and infiow conditions at the bottom face. 
Initially, molecular gas is static and has an ideal gas equation of state (7 =5/3), a number 
density of lOOcms"^ and a temperature of 10000 K. The magnetic field is helical, centrally 
localized and given by the vector potential (in cylindrical coordinates) 

A - / r) + l){cos{2 z) + 1)^ + f{cos{2 r) + l){cos{2 z) + 1)^, for r,z< re; 

"^^^'^^-1 0, forr,z>re, 

(1) 

where Tg = 30 AU and a =40 has units of length and determines the ratio of toroidal to 
poloidal magnetic fiuxes. The magnetic pressure exceeds the thermal pressure inside the 
magnetized region. 

Source terms continually inject magnetic or kinetic energy at cells r, z < rg. We carry 
out 4 simulations: an adiabatic, a cooling, a rotating PFD jet, and a hydro dynamical jet. 
The latter is constructed to have same time average propagation speed and energy fiux as 
the adiabatic PFD jet. The cooling tables of Dalgarno & McCray ()1972p were used for the 
cooling case. For the rotating case we allowed the gas and magnetic fields within r, z < rg to 
move in Keplerian rotation. 

3.2 Results 

Magnetic pressure pushes field 
lines and plasma upward, 
forming magnetic cavities with 
low density (Fig. 3, all but 
right panel). The adiabatic 
case is the most stable. PFD 
jets decelerate more quickly 
relative to our hydro jet. This 
results because the PFD and 
hydro jets have the same in- 
jected energy fiux, but the 
PFD case produces not only 
axial but radial expansion. 
The pre-collimated hydro can 
only expand via a much lower 
thermal pressure. Thus all of 
the energy fiux in the hyrdo- 
case for our set up is more ef- 
ficiently directed to axial me- 
chanical power. In principle. Figure 3: Logarithmic density maps of the adiabatic (1st col- 
our hydro case can emulate umn), rotating (2nd column) and cooling (3rd column) PFD 

jets. Hydrodynamic jet (4th column). From top to bottom 

the time is 42, 84 and 118 yr. 
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the asymptotic propagation regime of a jet that was magneto-centrifugahy launched (e.g. 
Blackman 2007) which is distinct from a PFD jet. 

The PFD jet cores are thin and unstable, whereas the hydro jet beam is thicker, smoother 
and stable. The PFD jets are sub-Alfvenic. Their cores are confined by magnetic hoop stress, 
while their surrounding cavities are collimated by external thermal pressure. PFD jets carry 
high axial currents which return along their outer contact discontinuity. 

We see that the inner regions of the PFD jets, just outside the core, are low beta plasma 
columns in which the axial magnetic field dominates over the toroidal one, \B(j)/ Bz\ <C 1. The 
columns' instability condition is given by 

^|>l(/3.-l)fcOe,|, (2) 

where j3z = 2iiqP/B1^ jhq is the magnetic vacuum permeability, P is the plasma thermal 
pressure and is the characteristic wavelength of the current-driven perturbations. We 
find that the cooling jet shows /^^ ~ 1 from early in the simulation (Fig. 4), and thus does not 
have sufficient thermal energy to damp the magnetic pressure kink perturbations. The latter 
grow exponentially. 




Adiabatic Rotating Cooling 

Figure 4: Plasma velocities along the axes of the PFD jets after expanding for 84 yr. The 
soHd, dashed, dotted, dot-dash and dot-dash thick lines represent the Alfven speed, the sound 
speed, Vz^Vx and Vy^ respectively. Each velocity unit represents 9.1kms~^. 

A different path to instability operates for the rotating (non-cooling) case. For this case, 
rotation at the base of the jet causes a slow ampliffcation of the toroidal magnetic ffeld. Hence 
the left hand side of equation (|2|) increases slowly, and so do the kink mode perturbations. 
The rotating jet is not completely destroyed by these perturbations, and their amplitude is 
about twice the radius of the central jet (Fig 5), in agreement with the Kruskal-Shafranov 
criterion (Kruskal et al. I1958( Shafranov ll958p . 

In Figure 4 we show profiles of the relevant velocities of the PFD jets along the jet 
axis, as a function of cooling and rotation. We also followed these in time. During the stable 
propagation phase, we find that the jets are mostly sub-Alfvenic and trans-sonic, independent 
of cooling or rotation. Fast-forward compressive MHD (FF) and transmitted hydrodynamic 
shocks are evident in the ambient medium, ahead of the jets. The FF shocks steepen in time. 
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whereas the hydrodynamic shocks are quickly dissipated in the coohng case. In contrast, the 
adiabatic and rotating cases show regions within the lower half of the jets, where the sound 
speed is super- Alfvenic. Such regions are bounded by the reverse and the forward slow- modes 
of compressive MHD waves, and characterized by high thermal to magnetic pressure ratios. 



4 Conclusions Magnetic field strength [/iG] 

PFD jets can be produced 15.0^5.00 35.00 45.00 5^00 

in pulsed power labora- 
tory facilities and help to 
understand the physics of 
astrophysical jets. The 
PFD "lab jets" are col- 
limated by hoop stress, 
and their outer magnetic 
cavities are collimated by 
external pressure. The 
axial magnetic component 
infiuences both the radial 
compression and the sta- 
bility of the jets. The 
outfiows evolve to become 
corrugated, but still col- 
limated, structures due 
to current-driven insta- 
bilities. Thin conduc- 
tion foil experiments pro- 
duce episodic jets and 
nested magnetic bubbles. 
Observations of multiple 
lobes in the radio galaxy 
B0925+420, may be con- 
sistent with such processes 
(Brocksopp et al. 2007). 

Our simulations show 
that PFD jet beams are 

lighter, slower and less stable than pre-collimated asymptotically hydrodynamic jets. In prac- 
tice the latter could represent the asymptotic propagation regimes of magneto-centrifugally 
launched jets, which are distinct from PFD in that PFD remain PFD out to much larger 
scales. We find that current-driven perturbations in PFD jets are amplified by both cooling 
and rotation for the regimes studied: Shocks and thermal pressure support are weakened by 
cooling, making the jets more susceptible to kinking. Rotation amplifies the toroidal magnetic 
field which also exacerbates the kink instability. Our simulations agree well with the models 
and experiments of Shibata & Uchida (jl986p and Lebedev et al. (|2005p , respectively. 
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Figure 5: Central magnetic field lines at t = 118yr. From left to 
right these are the adiabatic, the rotating and the cooling PFD 
jets, respectively. Bottom panels show a pole-on view. 
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